arXiv:1509.02372v3 [hep-ex] 29 Nov 2015 


EUROPEAN ORGANIZATION EOR NUCLEAR RESEARCH (CERN) 



CERN-PH-EP-2015-232 
LHCb-PAPER-2015-045 
September 8, 2015 


Forward production of T mesons in 
pp collisions at y/s = 7 and 8 TeV 


The LHCb collaboratioiiH] 


Abstract 

The production of T mesons in pp collisions at y/s = 7 and 8 TeV is studied with 
the LHCb detector using data samples corresponding to an integrated luminosity 
of 1 fb~^ and 2fb~^ respectively. The production cross-sections and ratios of 
cross-sections are measured as functions of the meson transverse momentum pT and 
rapidity y, for px < 30GeV/c and 2.0 < y < 4.5. 


Published in JHEP 1511(2015) 103 

© CERN on behalf of the LHCb collaboration, licence CC-BY-4.0 


^Authors are listed at the end of this paper. 




11 



1 Introduction 


In high energy hadron collisions, the prodnction of heavy qnarkoninm systems snch as 
the bb states (T(1S), T(2S) and T(3S), represented generically as T in the following) 
probes the dynamics of the colliding partons and provides insight into the non-perturbative 
regime of quantum chromodynamics (QCD). Despite many models that have been proposed, 
a complete description of heavy qnarkoninm production is still not available. 

The effective theory of non-relativistic QCD (NRQCD) provides the foundation 
for much of the current theoretical work. According to NRQCD, the production of heavy 
qnarkoninm factorises into two steps: a heavy quark-antiquark pair is hrst created at 
short distances, and subsequently evolves non-perturbatively into a qnarkoninm state. 
The NRQCD calculations include the colour-singlet (CS) and colour-octet (CO) matrix 
elements for the pertubative stage. The CS model [^|^, which provides a leading-order 
description of qnarkoninm production, underestimates the cross-section for single J/i|» pro¬ 
duction at the Tevatron at high px, where px is the component of the meson momentum 
transverse to the beam. To resolve this discrepancy, the CO mechanism was introduced |^. 
The corresponding matrix elements were determined from the high-px data, as the CO 
cross-section decreases more slowly with px than that predicted by the CS model. More re¬ 
cent higher-order calculations 00 show better agreement between CS predictions and 
the experimental data [^, reducing the need for large CO contributions. The production 
of T mesons in proton-proton (pp) collisions can occur either directly in parton scattering 
or via feed down from the decay of heavier bottomonium states, such as Xb 13 18 


or higher-mass T states, which complicates the theoretical description of bottomonium 


production 19,20 


The production of the T mesons has been studied using pp collision data taken 
at \/s = 2.76, 7 and 8 TeV by the LHCb 21-23, ALICE 24 , ATLAS 25 and CMS 26,2^ 


experiments in different kinematic regions. The existing LHCb measurements of these 
quantities were performed at y/s = 7TeV with a data sample collected in 2010 corre¬ 
sponding to an integrated luminosity of 25pb“^, and at y/s = 8TeV for early 2012 data 
using 50pb~^. Both measurements were differential in px and y of the T mesons in 
the ranges 2.0 < p < 4.5 and Px < 15 GeV/c. Based on these measurements, an increase 
of the production cross-section in excess of 30% between y/s = 7 and 8 TeV was ob¬ 
served, which is larger than the increase observed for other qnarkoninm states such as 
the J/ij; [^,28 


and larger than the expectations from NRQCD 11 


In this paper we report on the measurement of the inclusive production cross-sections of 
the T states at ^/s = 7 and 8 TeV and the ratios of these cross-sections. The T cross-section 
measurement is performed using a data sample corresponding to the complete LHCb Run 1 
data set with integrated luminosities of 1 fb“^ and 2 fb“^, accumulated at y/s = 7 and 8 TeV, 
respectively. These samples are independent from those used in the previous analyses |22p^ . 
The increased size of the data sample results in a better statistical precision and allows 
the measurements to be extended up to px values of 30GeV/c. 
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2 Detector and simulation 


The LHCb detector |29[|30| is a single-arm forward spectrometer covering the pseudorapidity 
range 2 < p < 5, designed for the study of particles containing b or c quarks. The detector 
includes a high-precision tracking system consisting of a silicon-strip vertex detector 
surrounding the pp interaction region, a large-area silicon-strip detector located upstream 
of a dipole magnet with a bending power of about 4 Tm, and three stations of silicon-strip 
detectors and straw drift tubes placed downstream of the magnet. The tracking system 
provides a measurement of momentum, p, of charged particles with a relative uncertainty 
that varies from 0.5% at low momentum to 1.0% at 200GeV/c. The minimum distance 
of a track to a primary vertex, the impact parameter, is measured with a resolution of 
(15 -|- 29 /pt) where px is in GeV/c. Different types of charged hadrons are distinguished 
using information from two ring-imaging Cherenkov detectors. Photons, electrons and 
hadrons are identified by a calorimeter system consisting of scintillating-pad and preshower 
detectors, an electromagnetic calorimeter and a hadronic calorimeter. Muons are identified 
by a system composed of alternating layers of iron and multiwire proportional chambers 


31 


The online event selection is performed by a trigger 32 , which consists of a hardware stage, 
based on information from the calorimeter and muon systems, followed by a software stage, 
which applies a full event reconstruction. At the hardware stage, events for this analysis 
are selected by requiring dimuon candidates with a product of their px values exceeding 
1.7 (2.6) (GeV/c)^ for data collected at a/s = 7 (8) TeV. In the subsequent software trigger, 
two well-reconstructed tracks are required to have hits in the muon system, px > 500 MeV/c, 
p > 6GeV/c and to form a common vertex. Only events with a dimuon candidate with 
a mass > 4.7 GeV/c^ are retained for further analysis. In the offline selection, trigger 

decisions are associated with reconstructed particles. Selection requirements can therefore 
be made on the trigger selection itself and on whether the decision was due to the signal 
candidate, the other particles produced in the pp collision, or a combination of both. 

In the simulation, pp collisions are generated using Pythia6 33 with a specific 
LHGb configuration |^. Decays of hadronic particles are described by EvtGen (^ , 
in which final-state radiation is generated using Photos j^. The interaction of the gen¬ 
erated particles with the detector, and its response, are implemented using the Geant4 
toolkit 


37 as described in Ref. 38 


3 Selection and cross-section determination 


The event selection is based on the criteria described in the previous LHCb T analy¬ 
ses 21 23 but slightly modified to improve the signal-to-background ratio. It includes 


selection criteria that ensure good quality track reconstruction 39 , muon identification 40 


and the requirement of a good fit quality for the dimuon vertex, where the associated 
primary vertex position is used as a constraint in the fit [^. In addition, the muon can¬ 
didates are required to have 1 < px < 25 GeV/c, 10 < p < 400 GeV/c and pseudorapidity 
within the region 2.0 < ri < 4.5. 
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The differential cross-section for the production of an T meson decaying into a muon 
pair is 


X 


d^ 


-ff(pp ^ TX) = 


-O' 




( 1 ) 


dpT dy ' Ap^Ay Ap^Ay C 

■where By is the branching fraction of the T —decay, Ay and Apx are the rapidity 
and pt bin sizes, is a production cross-section for T —events in the given 

(Pt,2/) bin, is the efficiency-corrected number of T —)■ p+p~ decays and C is 

the integrated luminosity. Given the sizeable uncertainty on the dimuon branching 


fractions of the T mesons 42 , the measurement of the production cross-section multiplied 


by the dimuon branching fraction is presented, as in previous LHCb measurements 21-23 


A large part of the theoretical and experimental uncertainties cancel in the ratios of 
production cross-sections of various T mesons, defined for a given (pt,p) bin as 


T(iS)^^+^- 
^bin 


N- 


T(iS)-s>n-l-n- 




(J 


bin 


N- 




( 2 ) 


The evolution of the production cross-sections as a function of pp collision energy is studied 
using the ratio 


ff 


bin 


^8/7 = 


^=8 TeV 


a 


bin 


( 3 ) 


^=7 TeV 

The signal yields A^T-s>^+^- in each (pt,p) bin are determined from an unbinned 
extended maximum likelihood fit to the dimuon mass spectrum of the selected candidates 
within the range 8.5 < m^+^- < 12.5 GeV/c^. The correction for efficiency is embedded 
in the ht procedure. Each dimuon candidate is given a weight calculated as 1/£*“*, where 
is the total efficiency, which is determined for each T —>■ p'^p” candidate as 


-.tot 


_ ^rec&sel ^ ^trg ^ ^|xID 


( 4 ) 


where is the reconstruction and selection efficiency, e*’’® is the trigger efficiency and 

is the efficiency of the muon identification criteria. The efficiencies and e*'"® are 

determined using simulation, and corrected using data-driven techniques to account for 
small differences in the muon reconstruction efficiency between data and simulation 39,40 


The efficiency of the muon identification criteria is measured directly from data using 
a large sample of low-background J/ij) —)■ p^p“ events. All efficiencies are evaluated as 
functions of the muon and dimuon kinematics. The mean total efficiency reaches 
a maximum of about 45% for the region 15 < px < 20 GeV/c, 3.0 < y < 3.5, and drops 
down to 10% at high px and large p, with the average efficiency being about 30%. 

In each (px, y) bin, the dimuon mass distribution is described by the sum of three 
Grystal Ball functions 43 , one for each of the T(1S), T(2S) and T(3S) signals, and 


the product of an exponential function with a second-order polynomial for the combinato¬ 
rial background. The mean value and the resolution of the Crystal Ball function describing 
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Figure 1: Efficiency-corrected dimuon mass distributions for (left) ^/s = 7TeV and 
(right) y/s = 8TeV samples in the region 3 < pT < 4GeV/c, 3.0 < y < 3.5. The thick dark 
yellow solid curves show the result of the hts, as described in the text. The three peaks, shown 
with thin magenta solid lines, correspond to the T(1S), T(2S) and T(3S) signals (left to right). 
The background component is indicated with a blue dashed line. To show the signal peaks clearly, 
the range of the dimuon mass shown is narrower than that used in the ht. 


the mass distribution of the T(1S) meson are free ht parameters. For the T(2S) and 
T(3S) mesons the mass differences m(T(2S)) — m(T(lS)) and m(T(3S)) — m(T(lS)) are 
hxed to the known values while the resolutions are hxed to the value of the reso¬ 
lution of the T(1S) signal, scaled by the ratio of the masses of the T(2S) and T(3S) to 
the T(1S) meson. The tail parameters of the Crystal Ball function describing the radiative 
tail are hxed from studies of simulated samples. 

The hts are performed independently on the efficiency-corrected dimuon mass dis¬ 
tributions in each {pT,y) bin. As an example, Fig. shows the results of the hts in 
the region 3 < pt < 4GeV/c and 3.0 < y < 3.5. For each bin the position and the res¬ 
olution of the T(1S) signal is found to be consistent between ^/s = 7 and 8TeV data 
sets. The resolution varies between 33 MeV/c^ in the region of low pt and small rapidity 
and 90MeV/c^ for the high pt and large y region, with the average value being close 
to 42MeV/c^. The total signal yields are obtained by summing the signal yields over all 
(Pt,p) bins and are summarised in Table 

4 Systematic uncertainties 

The systematic uncertainties are summarised in Table separately for the measurement 
of the cross-sections and of their ratios. 

The uncertainty related to the mass model describing the shape of the dimuon 
mass distribution is studied by varying the ht range and the signal and background 
parametrisation used in the ht model. The ht range is varied by moving the upper edge 
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Table 1: Efficiency-corrected signal yields for data samples accumulated at -y/s = 7 and 8TeV 
summed over the full kinematic range pT < 30GeV/c, 2.0 < y < 4.5. The uncertainties are 
statistical only. 


Vs = 7TeV Vs = 8TeV 

(2639.8 ± 3.7) ■ lO^ (6563.1 ± 6.3) ■ 10^ 

(667.3 ± 2.2) ■ lO^ (1674.3 ± 3.5) ■ 10^ 

A^t( 3S)^^+^- (328.8 ± 1.5) ■ lO^ (786.6 ± 2.6) ■ 10^ 

Table 2: Summary of relative systematic uncertainties (in %) for the differential production 
cross-sections, their ratios, integrated cross-sections and the ratios ^ 8 / 7 - The ranges indicate 
variations depending on the {pT,y) bin and the T state. 


Source 

^bin ^1,3 


'^8/7 

Fit model and range 

0.1-4.8 0.1-2.9 

0.1 

— 

Efficiency correction 

0.2-0.6 0.1-1.1 

0.4 

— 

Efficiency uncertainty 

0.2-0.3 — 

0.2 

0.3 

Muon identihcation 

0.3-0.5 — 

0.3 

0.2 

Data-simulation agreement 




Radiative tails 

1.0 — 

1.0 

— 

Selection efficiency 

1.0 0.5 

1.0 

0.5 

Tracking efficiency 

0.5 0(2x0.4) — 

0.5© (2 X 0.4) 

— 

Trigger efficiency 

2.0 — 

2.0 

1.0 

Luminosity 

1.7(y^= 7TeV) 

1.2(7^ = 8 TeV) 

1.7(y^ = 7TeV) 
1.2(7^ = 8 TeV) 

1.4 


from 12.5 to 11.5GeV/c^; the degree of the polynomial function used in the estimation of 
the background is varied between zeroth and the third order. Also the tail parameters of 
the Crystal Ball function are allowed to vary in the ht. In addition, the constraints on 
the difference in the T signal peak positions are removed for all bins with high signal yields. 
The maximum relative difference in the number of signal events is taken as a systematic 
uncertainty arising from the choice of the £t model. 

As an alternative to the determination of the signal yields from efficiency-corrected data, 
the method employed in Ref. 21 is used. In this method the efficiency-corrected yields 


for each (pt, 2/) bin are calculated using the sPlot technique [^. The difference between 
this method and the nominal one is taken as a systematic uncertainty on the efficiency 
correction. 

Reconstruction, selection and trigger efficiencies in Eq. (|^ are obtained using simulated 
samples. The uncertainties due to the hnite size of these samples are propagated to 
the measurement using a large number of pseudoexperiments. The same technique is used 
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for the propagation of the uncertainties on the muon identihcation efficiency determined 
from large low-background samples of J/i|) —)■ decays. 

Several systematic uncertainties are assigned to account for possible imperfections in 
the simulated samples. The possible mismodelling of the bremsstrahlung simulation for 
the radiative tail and its effect on the signal shape has been estimated in previous LHCb 
analyses 23 and leads to an additional uncertainty of 1.0% on the cross-section. 


Good agreement between data and simulation is observed for all variables used in 
the selection. The small differences seen would affect the efficiencies by less than 1.0%, 
which is conservatively taken as the systematic uncertainty to account for the disagreement 
between data and simulation. 

The efficiency is corrected using data-driven techniques to account for small differences 
in the tracking efficiency between data and simulation 39,40 . The uncertainty in the cor¬ 


rection factor is propagated to the cross-section measurement using pseudoexperiments 
and results in a global 0.5% systematic uncertainty plus an additional uncertainty of 
0.4% per track. 

The systematic uncertainty associated with the trigger requirements is assessed by 
studying the performance of the dimuon trigger for T(1S) events selected using the single 
muon high-px trigger 32 in data and simulation. The comparison is perfomed in bins of 


the T(1S) meson transverse momentum and rapidity and the largest observed difference of 
2 . 0 % is assigned as the systematic uncertainty associated with the imperfection of trigger 
simulation. 

The luminosity measurement was calibrated during dedicated data taking periods, 
using both van der Meer scans 45 and a beam-gas imaging method 46,47 . The abso¬ 


lute luminosity scale is determined with 1.7 (1.2)% uncertainty for the sample collected 
at -^/s = 7 ( 8 ) TeV, of which the beam-gas resolution, the spread of the measurements and 
the detector alignment are the largest contributions 47 49 . The ratio of absolute lumi¬ 


nosities for samples accumulated at ^/s = 7 and 8 TeV is known with a 1.4% uncertainty. 

The total systematic uncertainty in each {pT,y) bin is the sum in quadrature of 
the individual components described above. For the integrated production cross-section 
the systematic uncertainty is estimated by taking into account bin-to-bin correlations. 
Several systematic uncertainties cancel or signihcantly reduce in the measurement of 
the ratios and ^ 8/71 as shown in Table 

The production cross-sections are measured at centre-of-mass energies of 7 and 8 TeV, 


where the actual beam energy for pp collisions is known with a precision of 0.65% 50 


Assuming a linear dependence of the production cross-section on the pp collision energy, 
and using the measured production cross-sections at a/s = 7 ( 8 ) TeV, the change in the pro¬ 
duction cross-section due to the imprecise knowledge of the beam energy is estimated 
to be 1.4 (1.2)%. The effect is strongly correlated between a/s = 7 and 8 TeV data and 
will therefore mostly cancel in the measurement of the ratio of cross-sections at the two 
energies. 

The efficiency is dependent on the polarisation of the T mesons. The polarisation of 
the T mesons produced in pp collisions at ^/s = 7 TeV at high px and central rapidity has 
been studied by the CMS collaboration in the centre-of-mass helicity, Collins-Soper 52 


6 
















and the perpendicular helicity frames. No evidence of significant transverse or longitudinal 
polarisation has been observed for the region 10 < pt < 50GeV/c, \y\ < 1.2. Therefore, 
results are quoted under the assumption of unpolarised production of T mesons and no cor¬ 
responding systematic uncertainty is assigned on the cross-section. Under the assumption 
of transversely polarised T mesons with = 0.2 in the LHCb kinematic region]^ the total 
production cross-section would result in an increase of 3%, with the largest local increase 
of around 6% occuring in the low region (pt < 3 GeV/c), both for small {y < 2.5) and 
large (p > 4.0) rapidities. 


5 Results 


The double-differential production cross-sections multiplied by the dimuon branching 


fractions for the T mesons are shown in Fig. T 
in (pt,p) bins are presented in Tables 


CT, 




bin 


ae corresponding production cross-section 
3, 4 and for ^/s = 7 TeV and Tables 6 
and 1^ for -y/s = 8 TeV. The cross-sections integrated over p as a function of pt anc 
integrated over px as a function of rapidity are shown in Figs. |^and|^ respectively. 

The transverse momentum spectra are £t using a Tsallis function (5^ 


da 


PxdpT 


oc 1 -h 


rpkin 

nT 


(5) 


where = 


= 


X + Px ~ is the transverse kinetic energy, the power n and the tem¬ 
perature parameter T are free parameters, and mr is the known mass of a T meson (4^ . 
This function has a power-law asymptotic behaviour oc p^" for high px as expected for 
hard scattering processes. It has been successfully applied to ht px spectra 5^57 in wide 
ranges of particle species, processes and kinematics. A ht with the Tsallis distribution 
for the range 6 < px < 30GeV/c is superimposed on the differential cross-sections in 
Fig-i The ht quality is good for all cases. The htted values of the parameters n and 
T are listed in Table The parameter n for all cases is close to 8, compatible with 
the high px asymptotic behaviour expected by the CS model |3t|H [58| - |60] . The temperature 
parameters T show little dependence on y/s and increase with the mass of T state. 

The shapes of the rapidity spectra are compared with the GO model prediction in 
the region 2.5 < y < 4.0 and are htted using the function given by Eq. (1) of Ref. (^ , 
with free normalisation constants. The ht result, as well as the extrapolation to the full 
kinematic range 2.0 < y < 4.5, is presented in Fig. The quality of the ht is good for 
all cases. 

The integrated production cross-sections multiplied by the dimuon branching fractions 
in the full range px < 30GeV/c and 2.0 <y < 4.5 at ^/s = 7 and 8 TeV are reported 


in Table IT where the hrst uncertainties are statistical and the second systematic. The same 
measurements are also shown integrated over the reduced range px < 15 GeV/c in the same 
rapidity range, to allow the comparison with previous measurements 1^1^. 


^The CMS measurements for T(1S) mesons are consistent with small transverse polarisation in 
the helicity frame with the central values for the polarisation parameter 0 < A ,5 < 0.2 51 . 
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(bottom) T(3S) at (left) yG = 7TeV and (right) yG = 8 TeV. The error bars indicate the sum 
in quadrature of the statistical and systematic uncertainties. The rapidity ranges 2.0 < y < 2.5, 
2.5 < y < 3.0, 3.0 < y < 3.5, 3.5 < y < 4.0 and 4.0 < y < 4.5 are shown with red filled circles. 


blue open squares, cyan downward triangles, magenta upward triangles and green diamonds, 
respectively. Some data points are displaced from the bin centres to improve visibility. 
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Table 3: Production cross-section ^ [pb] in {pT,y) bins for ^/s = 7TeV. The first uncertainties are statistical and 

the second are the uncorrelated component of the systematic uncertainties. The overall correlated systematic uncertainty is 3.1% 
and is not included in the numbers in the table. The horizontal lines indicate bin boundaries. 
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Table 4: Production cross-section ^ [pb] in {pT,y) bins for y/s = 7TeV. The first uncertainties are statistical and 

the second are the uncorrelated component of the systematic uncertainties. The overall correlated systematic uncertainty is 3.1% 
and is not included in the numbers in the table. The horizontal lines indicate bin boundaries. 
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Table 5: Production cross-section ^ [pb] in {pT,y) bins for y/s = 7TeV. The first uncertainties are statistical and 

the second are the uncorrelated component of the systematic uncertainties. The overall correlated systematic uncertainty is 3.1% 
and is not included in the numbers in the table. The horizontal lines indicate bin boundaries. 
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Table 6: Production cross-section ^ [pb] in {pT,y) bins for y/s = 8TeV. The first uncertainties are statistical and 

the second are the uncorrelated component of the systematic uncertainties. The overall correlated systematic uncertainty is 2.8% 
and is not included in the numbers in the table. The horizontal lines indicate bin boundaries. 
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Table 7: Production cross-section ^ [pb] in {pT,y) bins for y/s = 8TeV. The first uncertainties are statistical and 

the second are the uncorrelated component of the systematic uncertainties. The overall correlated systematic uncertainty is 2.8% 
and is not included in the numbers in the table. The horizontal lines indicate bin boundaries. 
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Table 8: Production cross-section ^ [pb] in {pT,y) bins for y/s = 8TeV. The first uncertainties are statistical and 

the second are the uncorrelated component of the systematic uncertainty. The overall correlated systematic uncertainty is 2.8% and 
is not included in the numbers in the table. The horizontal lines indicate the bin boundaries. 
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Figure 3: Differential cross-sections in the range 2.0 < y < 4.5 for (red solid 

circles) T(1S), (blue open squares) T(2S) and (green solid diamonds) T(3S) mesons for 
(left) y/s = 7TeV and (right) y/s = 8TeV data. The curves show the fit results with the Tsal- 
lis function in the range 6 < px < 30 GeV/c. The data points are positioned in the bins according 
to Eq. (6) in Ref. 


61 



Figure 4: Differential cross-sections in the range px < 30GeV/c for (red solid 

circles) T(1S), (blue open squares) T(2S) and (green solid diamonds) T(3S) mesons for 
(left) y/s = 7 TeV and (right) y/s = 8 TeV data. Thick lines show ht results with the CO model 
predictions from Ref. in the region 2.5 < y < 4.0, and dashed lines show the extrapolation 

to the full region 2.0 < y < 4.5. The data points are positioned in the bins according to Eq. (6) 
in Ref. 
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The ratios of integrated production cross-section <^^ 8/7 nre presented in Table [TT] for 
the full (px < 30 GeV/c) and reduced (px < 15 GeV/c) ranges. The results for the reduced 
range are consistent with the previous measurements, confirming the increase of the bot- 
tomonium production cross-section of approximately 30% when the centre-of-mass energy 
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Table 9: Results of the fits to the transverse momentum spectra of T mesons using the Tsallis 
function in the reduced range 6 < px < 30GeV/c. 




T [GeV] 

n 

T(1S) 

7 TeV 

1.19 ±0.04 

8.01 ±0.33 

8 TeV 

1.20 ±0.04 

7.71 ±0.27 

T(2S) 

7 TeV 

1.33 ±0.05 

7.57 ±0.41 

8 TeV 

1.37 ±0.05 

7.53 ± 0.34 

T(3S) 

7 TeV 

1.53 ±0.07 

7.85 ±0.56 

8 TeV 

1.63 ±0.06 

8.23 ±0.51 


Table 10: The production cross-section (in pb) for T mesons in the full kinematic 

range px < 30GeV/c(left two columns), and reduced range px < 15 GeV/c (right two columns), 
for 2.0 < y < 4.5. The first uncertainties are statistical and the second systematic. 


Px < 30 GeV/c 


Px < 15 GeV/c 


s = 7TeV 


s = 8 TeV 


s = 7TeV 


s = 8 TeV 



2510±3±80 

3280 ± 3 ± 100 

2460 ± 3 ± 80 

3210 ± 3 ± 90 


635 ± 2 ± 20 

837±2±25 

614 ± 2 ± 20 

807 ± 2 ± 24 


313±2± 10 

393 ± 1 ± 12 

298 ± 1 ± 10 

373 ± 1 ± 11 

increases from 

00 

o 

22 23 . 




The ratios <^^ 8/7 as a function of px integrated over the region 2.0 < y < 4.5 are 
shown in Fig. |^. The ratios are htted with a linear function. The ht quality is good, 
with a p-value exceeding 35% for all cases, and the slopes are found to be 10.8 ± 0.6, 
9.5 ± 1.2 and 9.8 ± 1.6 (in units of 10“^/ (GeV/c)) for T(1S), T(2S) and T(3S), respec¬ 
tively. The measurements are compared with the NRQCD theory predictions in 
the same kinematic range, where only uncertainties from the CO long distance matrix 
elements are considered since most other uncertainties are expected to cancel in the ratio. 
The theory predictions are independent on the T state and are consistently lower than 
the measurements. 

The ratio ^ 8/7 as a function of rapidity, integrated over the region px < 30 GeV/c is 
shown in Fig. [^. The ratios are compared with the expectations from the CO mecha¬ 
nism 62 ,^ with normalisation factors hxed from the hts of Fig. The trend observed 
in data does not agree with the pure CO model. It can be noted that also for open 
beauty hadrons the differential cross-sections exhibit a larger rise as a function of a/s 


at smaller rapidities 54 , while the FONLL calculations 64 predict this behaviour towards 
larger rapidity. 

The ratios at a/s = 7 and 8 TeV are reported in Fig. and Tables and 
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Table 11: The ratio of production cross-sections for T mesons at ^/s = 8 to that at ^/s = 7TeV 
in the full kinematic range px < 30 GeV/c (left) and reduced range px < 15 GeV/c (right) 
for 2.0 < y < 4.5. The first uncertainties are statistical and the second systematic. 


Pt < 30 GeV/c 


Pt < 15 GeV/c 


T(1S) 1.307 ±0.002 ±0.025 

T(2S) 1.319 ±0.005 ±0.025 

T(3S) 1.258 ±0.007 ±0.024 


1.304 ±0.002 ±0.024 
1.315 ±0.005 ±0.024 
1.254 ±0.007 ±0.023 



3 

l> 

.^00 



Figure 5: Ratios of the differential cross-sections (left) and (right) 

at ^/s = 8 and 7TeV for (red solid circles) T(1S), (blue open squares) T(2S) and (green solid 
diamonds) T(3S). On the left hand plot, the results of the fit with a linear function are shown with 
straight thin red solid, blue dotted and green dashed lines. In the same plot, the next-to-leading 
order NRQCD theory predictions 11 are shown as a thick line. On the right hand plot, the curved 

with 


red solid, blue dotted and greed dashed lines show the CO model predictions (62,63 


the normalisation hxed from the fits in Fig. |^for T(1S), T(2S) and T(3S) mesons, respectively. 
Some data points are displaced from the bin centres to improve visibility. 


as a function of px for different rapidity bins. The same ratios as a function of px 
integrated over rapidity, and as a function of y integrated over px, are shown in Fig. 


The ratios show little dependence on rapidity and increase as a function of px, 

and CMS (261 


ij 


ATLAS 
at 


25 


s = 7 and 8TeV are 


in agreement with previons observations by LHCb 22,23 
at a/s = 7TeV. The ratios of integrated cross-sections ^ 
reported in Table 16, for the full and the rednced px kinematic regions. All ratios 
agree with previons LHCb measurements. The ratio agrees with the estimates of 
0.27 from Refs. 63,65 , while signihcantly exceeds the expected value of 0.04 63 65 


but agrees with the range 0.14 — 0.22, expected for the hypothesis of a large admixture of 


a hybrid quarkonium state in the T(3S) meson state 65 
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Figure 6: The production ratios for (top) T(2S) to T(1S), (middle) T(3S) to T(1S), 
and (bottom) T(3S) to T(2S), measured with data collected at (left) \/s = 7TeV and 
(right) y/s = 8 TeV. The error bars indicate the sum in quadrature of the statistical and system¬ 
atic uncertainties. The rapidity ranges 2.0 < y < 2.5, 2.5 < y < 3.0, 3.0 < y < 3.5, 3.5 < y < 4.0 
and 4.0 < y < 4.5 are shown with red circles, blue squares, cyan downward triangles, magenta 
upward triangles and green diamonds, respectively. Some data points are displaced from the bin 
centres to improve visibility. 
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Table 12: The ratio ^ 2,1 for y/s = 7TeV. The first uncertainties are statistical and the second are the uncorrelated component of 
the systematic uncertainties. The overall correlated systematic uncertainty is 0.7% and is not included in the numbers in the table. 
The horizontal lines indicate bin boundaries. 
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Table 13: The ratio ^34 for y/s = 7TeV. The first uncertainties are statistical and the second are the uncorrelated component of 
the systematic uncertainties. The overall correlated systematic uncertainty is 0.7% and is not included in the numbers in the table. 
The horizontal lines indicate bin boundaries. 
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Table 14: The ratio ^ 2,1 for y/s = 8 TeV. The first uncertainties are statistical and the second are the uncorrelated component of 
the systematic uncertainties. The overall correlated systematic uncertainty is 0.7% and is not included in the numbers in the table. 
The horizontal lines indicate bin boundaries. 
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Table 15: The ratio ^34 for y/s = 8 TeV. The first uncertainties are statistical and the second are the uncorrelated component of 
the systematic uncertainties. The overall correlated systematic uncertainty is 0.7% and is not included in the numbers in the table. 
The horizontal lines indicate bin boundaries. 
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Figure 7: The production ratios (red solid circles) ^ 2 , 1 , (blue open squares) and (green 
solid diamonds) ^ 3^2 for (left) = 7TeV and (right) -y/s = 8TeV data, integrated over 
the (top) 2.0 < y < 4.5 region and (bottom) pT < 30GeV/c region. 


6 Summary 


The forward production of T mesons is studied in pp collisions at centre-of-mass energies 
of 7 and 8 TeV using data samples corresponding to integrated luminosities of 1 fb~^ and 
2 fb~^ respectively, collected with the LHCb detector. The double differential production 
cross-sections are measured as a function of meson transverse momenta and rapidity 
for the range < 30GeV/c, 2.0 < y < 4.5. The measured increase in the production 
cross-sections of T mesons between ^/s = 8 and 7 TeV significantly exceeds theory expecta¬ 
tions and confirms the previous LHCb observations (^ ,23 . For the region px < 15GeV/c 
the results agree with the previous measurements [22,23 , and supersede them. 
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Table 16: The ratios in the full kinematic range px < 30GeV/c and in the reduced 
range px < 15GeV/c for 2.0 < y < 4.5. The first uncertainties are statistical and the second 
systematic. 


v/s = 7TeV v^ = 8TeV 


Px < 30 GeV/c 


^2,1 0.253 ± 0.001 ± 0.004 0.255 ± 0.001 ± 0.004 

^34 0.125 ±0.001 ±0.002 0.120 ±0.000 ±0.002 

0.493 ± 0.003 ± 0.007 0.470 ± 0.002 ± 0.007 


px < 15 GeV/c 


^2,1 0.249 ± 0.001 ± 0.004 0.251 ± 0.001 ± 0.004 

^34 0.121 ±0.001 ±0.002 0.116 ±0.000 ±0.002 

M 3,2 0.485 ± 0.003 ± 0.007 0.463 ± 0.002 ± 0.007 
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